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ON THE FLOW OF A COMPRESSIBLE FLUID BY THE HODOGRAPH METHOD 
I— UNIFICATION AND EXTENSION OF PRESENT-DAY RESULTS 

By I. E. Gabkick and Cabi. Kaplan 



SUMMARY 

Elementary bamc sohiions of the eguaiions oj motion oj a 
compressible fluid in the hodograph variables are developed and 
used to provide a basis for comparison, in the form of velocity 
correction formulas, of corresponding compressible and incom- 
pressible flows. The hnoum approximate results of Ohaplygin, 
von Kdrmdn and Tsien, Temple and Yarwood, and Prandtl 
and Olauerl are unified by means of the analysis of the 
present paper. Two new types of approximations, obtained 
from the basic solutions, are introduced; they possess certain 
desirable features of the other approximxiMons and appear 
preferable as a basis for extrapolation into the range of high 
stream Mach numbers and large disturbances to the main 
stream. Tables and flares giving velocity and pressure- 
coefficient correction factors are included in order to facilitate 
fhf practical application of the results. 

INTRODUCTION 

The present paper is concerned with a theoretical study 
of the hydrodynamical equations of a perfect compressible 
fluid in two dimensions, in which the so-called hodograph 
variables are used as the independent variables. It is hoped 
to achieve herein a unification of the present-day results 
obtained in this field and also to provide a working basis for 
further developments. The earliest contributors to the 
hodograph method for treating compressible fluids were 
Molenbroek (reference 1) and Chaplygin (reference 2). The 
remarkable work of Chaplygin on gas jets appeared in 
Russian in 1904 but remained relatively imnoticed. In 
recent years contributions to tbe hodograph method have 
been made chiefly by Demtchenko (reference 3), von X&rm&n 
(reference 4), Tsien (reference 6), Ringleb (reference 6), and 
Temple and Yarwood (reference 7). 

The chief reason, and perhaps the only reason, for pre- 
ferring the hodograph variables to the physical plane co- 
ordinates is that the equations of motion in the hodograph 
variables are linear. This simplification is achieved, how- 
ever, at the cost of more difficult boimdary conditions and 
at a loss of physical insight. The great simplification in the 
mathematics due to linearity nevertheless makes it desirable 



to pursue this line of attack as long as it appears profitable 
to do so. 

The mathematics for handling the flow equations re- 
ceived a substantial impetus by the work of Bers and 
Gelbart (reference 8), who developed a new function theory 
analogous to ordinary analytic function theory. The 
present paper utilizes the methods of this new function 
theory to develop certain functions essential to the compres- 
sible-flow problem. It is of historical interest that ideas 
similar to those of Bers and Gelbart were explored by the 
renowned mathematician Hilbert (reference 9) in the early 
part of this century but do not appear to have been further 
developed at the time. 

The material to be treated is conveniently separated into 
two parts. In part I, the present paper, basic particular 
solutions of the hodograph flow equations are developed and 
employed in imifying and ejctending the results obtained by 
Chaplygia, von K&rmdn, and Temple and Yarwood. The 
results obtained in part I are of immediate practical appHca- 
tion and are given in the form of tables and graphs of velocity 
and pressxu-e-coefficient correction factors. In part EE, 
general particular solutions of the hodograph flow equations 
are developed and discussed. The material in part 11, it is 
hoped, will lead to a method for handling the actual boundary 
problem of the flow of a compressible fluid past a prescribed 
body. 

ANALYSIS 

FLOW EQUATIONS OF AN INCOMPBBSSIBLE FLUID 

It is well known that the relations between the velocity 
potential <j> and the stream function ^ for the steady irrota- 
tional two-dimensional motion of a perfect incompressible 
fluid are 

^y 



(1) 



These equations are the Cauchy-Riemann equations and 
therefore 4>-\-i^ is an analytic function /(z) of the comples^ 
variable z=x-\-iy. 
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aer-i' 



The complex velocity or reflected velocity vector u—iv is 
obtained from the complex potential /(z) by differentiation. 
Thus, 

df(z) 



u—tv- 



26 

g-fiS+t log q) 



(2) 



where g is the magnitude of the velocity vector and 0 is the 
angle the vector makes with the positive direction of the 
a>axis. 

The variables 0 and q are sometimes referred to as "the 
hodograph variables." The flow equations in the variables 
0 and g can be readily derived by introducing d+i log g as the 
independent complex variable in place of x+iy. Then, in 
analogy with equation (1), 



dfl d log 2 



d0 



or 



& log 2" 

dip 

d2 2 



(3) 



(4) 



(5) 



These equations are known as the hodograph equations for 
the flow of an incompressible flmd. 

FLOW EQUATIONS OF A COMPBESSIBLE FLUID 

The equations corresponding to equation (1) are, for a 
compressible fluid, 

i>z~p Z)y 

b<t> po ^ 

by~ p dxj 

where p is the density of the fluid at any point (x,y) and po is 
a constant density, which for convenience is referred to a 
stagnation point. 

A short way to derive the hodograph equations for a 
compressible fluid, attributed to Molenbroek, is as follows: 

According to equations (5), with ''i="5~ ojid. 

d<l>+i — d4f= (u dx+v dy) +i (— » dx-\-u dy) 
p 



or 



= {31— iv) {dx-\-i dy) 
=qer^ dz 

dz=h''(d4-\-i^>P^ 



(6) 



It follows fronjL equation (6), by considering 6 and 2 ^ 
independent variables, that 



50" 2^ p dey 



and 



d2 2 V^S P 2>2/ 



Then, by assuming that p is a function of only 2 (equivalent 
to assuming that the pressure is a function of only tlie 
density), 

d2&«? L 2* dq -deya \d2d9^^ p i>2W 



and 



&e&2~<z v^a p ^2/ 2 vaoaa^ p a0&2/ 

Since, by continuity, these two expressions are identical, it 
follows that 

Hence, by equating real and imaginary parts, 

" ' be p bq 



d2~2 dq bB 



(7) 



These are the hodograph equations, first obtained by 
Molenbroek, for the flow of a compressible Quid and are 
independent of the form of the pressure-density relation. 
It is observed that, when p=po=Constant, equations (7) 
reduce to equations (4). Equations (7), in contrast with 
equations (5), are linear in the dependent variables. 

BEBNOUUJ'S EQUATION AND EQUATION OF STATE 

In the present section there is listed a collection of for- 
mulas and definitions necessary in the analysb. 
Bernoulli's equation for a compressible fluid is 



*''"+|2*=0 



(8) 



jpt P 
where 

p static pressure in fluid 

fo static pressure at stagnation point (2=0) 

p density of fluid 

2 magnitude of velocity of fluid 

The adiabatic relation between the pressure and the 
density is 

s.=.(£.y (9) 

Po \Po/ 

where 

y adiabatic index (approx. 1.4 for air) 
Po density of fluid at stagnation point (2=0) 
The local velocity of sound a is obtained from 



dp 



For the adiabatic case. 



2 V 
p 



(10) 
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From Bornoulli's equation (8) and from equations (9) 
and (10), the following relations may be obtained: 



1 

-if-i 



p=po[i-|(t-i)^ 

.=.o[l-|(.-l)|]^, 



(11) 



where Oo is the velocity of sound at stagnation point (2=0). 

From equations (11) for 7]>1, a maximmn velocity 
2=g„ is obtained for the limiting conditions p=p=a=0. 
Thus, 

=2|8ao* (12) 

where 

The fimdamental nondimensional speed variable, in 
general, is g/oo but it is found useful in the analysis to employ 
a nondimensional speed variable r defined as 



(13) 



For 'y>l, the range of the variable t is O^t^ 1. The value 
T=0 has a dual meaning; t=0 in the case of a compressible 
fluid corresponds to a stagnation point (2=0), or t=0 may 
mean tlie limiting case of an incompressible fluid (ao=oo). 

With the definitions of r and jS, equations (1 1) become 

a=ao(l-T)'/» 

p=Po(l — r)^ 



P=Po(l-r)«" 



(14) 



The local Mach number M=-f may be expressed in terms 



of the speed variable r in the following way: 

2j3r 
^1— T 

or, by solving for t in terms of M, 

W 



(15) 



(16) 



The value of t for which the local velocity of the fluid 
equals the local velocity of sound (Af=l) is given by 



2(3+1 



(17) 



In the case of imiform flow past a fixed bomdary, the 
pressure coefficient is defined as 

843110—30 19 



2 "'21 

where the subscript 1 refers to the undisturbed stream. The 
pressure coefficient for the incompressible case (M=0) is 



The pressure coefficient for the compressible case is 



1+ i+^(-y-i)M 



For 2=2, (so»tc), 



,C ^ 2 p+(7-l)Afia -^ 



For 2=2m (focttttm), 



iPvM^m 



(18a) 



(18b) 



(18c) 



(18d) 



BASIC SOLimONS OF HODOGEAPH EQUATIONS 



Consider the incompressible case represented by equa- 
tions (3) or (4). It is dear that <^=fl and ^=log 2 satisfy 
these equations. In fact, any convergent power series in 
to=5+i log 2 represents an analjiiic function of which the 
real and imaginary parts satisfy equations (3) or (4). The 
class of analytic functions in w (and the concept of analytic 
continuation) then yields all the particular solutions of these 
equations. 

The particular solution w=0+i log 2 can be obtained by 
means of an integration that is instructive in the generaliza- 
tion to the compressible case. It is well known that 

i?'(w)=J/(w)rfti; 

can be represented as the sum of two line integrals 

F(?B)=^ (Pde-Qdlogq)+ij (Q de+P rflog q) 

whore 

Thus, given a pair of functions P and Q that satisfy equations 
(3) or (4), this process yields another pair of solutions, 
namely, the real and the imagmary parts of Fiw). For 
example, if P=l and (2=0> 



, F{w) =w=6+i log q 
Again, if P=0 and ^=1, 

F(w) =iv}= —log q,+i6 



(19) 



(20) 



The physical interpretation of equations (19) and (20), 
considered as flow patterns, is of some interest in connection 
with later developments. It is clear that equations (19) 
and (20) represent a vortex and a source located at the 
origin, respectively. 
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The generalization to the compressible case of the fore- 
going elementary results was accomplished by Bers and 
Gelbart (reference 8) by means of simple yet fertile ideas. 
Bers and Gelbart treat equations of the form 



(21) 



and show as is read£y verified that, if P and Q are a pair 
of solutions, the real and imaginaiy parts of the following 
sum of line integrals 

J[P dd-i^m di+ij de+^ P da] (22) 

are also solutions of equations (21). • 

In particular, corresponding to the pair of solutions P=l 
and Q=Q, there is obtained 

and, for P=0 and Q=l, 

iW=i[0+ij-S2(a) dg\ (24) 

By repeated application of the process of integration, indi- 
cated by expression (22), a general set of particular solutions 
of equations (21) may be obtained. These particular solu- 
tions are discussed in part U; in the present paper, only the 
solutions given by equations (23) and (24) are needed. 

The general hodograph equations (7) are of the form of 
equations (21) with 



and 



Xi{2) = 



For the rest of this paper, the adiabatic pressure-density 
relation (9) is used. By means of equations (9) and (14) 
and the relation 

dg q 

obtained from the differential form of Bemoxilli's equa- 
tion (8), it follows that 



Xi(2) = 



o--ry 



(26) 



The evaluation of the integrals in equations (23) and (24) 
is mode unique by requiring that the results reduce to the 
incompressible case when the speed of sound is infinite 
(that is, when t=0. Then, 



where 
and 

where 



=log 2-h/(r) 

=log 2+ff(T) 
1 r'-r i-(2j3+l )T ndr 



(26) 



(27) 



and it is observed that the functions /(r) and gir) vanish 
for T=0. 

Equations (23) and (24) can be written in the form 
W=e+iL 

and 

iW=i(B+iL) 

It is important to note that, in the incompressible case, 
W and iW reduce to w and iw, since L and L reduce to log g. 
Thus, there are in the compressible case two basic functions 
L and L corresponding to the one function log 2 in the in- 
compressible case. It is of interest to mention that the 
fimctions W and iW, considered as flow patterns in a com- 
pressible fluid, can again be interpreted as a vortex and a 
source. 

EVALUATION OF FUNCTIONS /(t) AND j(t) FOK VARIODS VALUES OF P 

In general, the integrals in equations (26) and (27) repre- 
senting the functions /(t) and g{r) are expressible by infinite 
series. For the important case of air, however, with the 
adiabatic index 7 put equal to 1.4 instead of the usual value 
1.408, these functions can be obtained in closed forms. 
Thus, with /3=2.6, 

/(r)=|j;[(l-r)--l]^^ 

= 1 H-ry+l (l-r)'/* 

+(i-.)--f|-iogl±iiF^ 



1/2 



(28) 



and 



1 r-r i-6r ndr 



(1-t)»/«T3(1-t)''/» 



(1-r) 



1 , i+a-ry* 

iS-g-log 2 



(29) 



Table 1 contains values of /(t) and gir), and figiure 1(a) 
shows these functions plotted against r. Observe that/(T) 
and gir) ore weU-behaved functions in the range 0^t<1. 
In figure 1(b), these functions are plotted against the local 
Mach number M in the practical Speed range. 
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Other interesting cases for which the functions /(t) and 

3 

gir) can be expressed in closed forms are y=^, y=2, 7=2' 
and 7=— 1. For 7=00 (/3=0, o=a>, incompressible case), 



For 7=2 (j3=l), 



f(T)=ff(r)=0 
/(r) = -ir 



?«=l-lZr^-ilog (1-r) 



For 7=1 (|8=2) 



At) r+^T* 



,1 1 1+r 
ff(r)=o-o7T— 



^ log (1-r) 



'2 2 (l-ry 2 

For7=-r(/3=-i), 

y«=,W=_logi±i^f^^'^ 

For tbe isothermal case 7=l(/3=oo), the velocity of sound 
o=afl=Constant;and tihe functions / and g are obtained as 
infinite series in the ratio g/oo. Thus, in the limit p->o>, 

_1 r^°«/-^g ,] d{gja,) 
~2jo 7 



-2-1^ A.' 2«+»W7i! 



n-l 



and 



1 C^'^ 



<^(g/ao) 
2/ao 



= 2 Jo K'-y^ ~^J^ 



For arbitrary^values of 7 (or j8) the expressions iorfir) and 
(t), obtained with the aid of the binomial expansion, are 



/(.)4|;(-i)-g)^ 

= -| /3r+| 13(13- 1)t*- 
1^.1 



4 00= 



+3k(2-T)(;g)--.. 



and 



= |/3r-|/3(^+l)r»- • . . 

~ 4 ^ 32 ^ W/ 

The significant feature of this general result is that, if powers 
of qjoo higher than the third are neglected, 



(30) 



and does not involve explicitly the adiabatic index 7. This 
circumstance imderlies the present-day approximate methods 
for obtaining velocity and pressure-coefficient correction 
factors; in the following sections, this point is brought out 
more clearly. 

APPUCATION OF BASIC FUNCTIONS L AND £ 

In this section, the basic functions L and L are employed 
to set up relations between velocities in "con-esponding" 
compressible and incompressible flows. These relations arc 
of the nature of "stretching factors" or velocity correction 
formulas and contain the results of Chaplygin, von Kdrmdn, 
Temple and Tarwood, and Glauert and Prandtl. It is 
important to recognize at the outset that no single velocity 
correction formula can represent in an exact way the cor- 
respondence of flow patterns past a prescribed body in 
a compressible and an incompressible fluid. A single velocity 
correction formula is actually feasible in only two cases: 

(1) The stream Mach niunber is small (even though tho 
disturbance to the main stream due to tho presence of the 
body may be large) so that the compressible-flow pattern 
differs only slightly from the incompressible-flow pattern or 

(2) the disturbance to the main stream is vanishingly small 
(even though the stream Mach number may be high) so that 
the effect of the shape of the soMd boimdary is small. Tho 
various velocity correction formulas discussed in the present 
paper differ essentially only in the degree to which tho 
requirements of these two cases are satisfied. Despite their 
limitations, single velocity correction formvdas are extra- 
polated, in view of the lack of more rigorous solutions, into 
the range of large disturbances to the main stream and high 
Mach numbers. This extrapolation can be justified by 
further theoretical investigations and by comparison ■svith 
experimental resiflts. 

Consider again the corresponding pairs of functions 

w=6+i log 2 
W=e+iL 



(31) 



and 



i'w=i{d+i log 2) 



(32) 
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It laos previously been noted that the pairs of fiinctions in 
equations (31) and (32) denote respectively a vortex and a 
source in an incompressible and a compressible fluid. Each 
pair of functions can be employed to define a correspondence 
of flow patterns in which corresponding points are identified 
by the same values i<f>,i'). Thus, in the case of the vortex 
(equations (31)), 

^f'^=^^.=log g,=i 

where the subscripts i and c refer to the incompressible and 
to the compressible case, respectively. It follows that 

=^yw (33) 
Similarly, in the case of the source (equations (32)), 
*f=<^<.=— Iog2*=— Z 



and 



(34) 



At the end of the preceding section it was pointed out that, 
to a first approximation, the functions/(T) and g(,T) are equal. 
This fact implies that, to a first approximation, a single veloc- 
ity correction formxila is feasible. The assumption is now 
made that either equation (33) or equation (34) can be adopted 
to provide a correspondence of flow patterns in the case of 
unifoim flow past a body in an incompressible and a com- 
pressible fluid. With the xmdisturbed streams as con- 
venient references, the following nondimensional forms of 
equations (33) and (34) can be written: 



/(t) 

7^ 



and 



(36) 
(36) 



where the subscript 1 refers to the imdisturbed stream. 
The use of the undisturbed stream as reference in the non- 
dimensional form of the velocity correction formula was 
introduced by Tsien in reference 5, where also the details of 
the von Kdrmdn approximation are developed. It is shown 
in the following section that either of equations (35) or (36) 
contains the result of Chaplygin, von Edrmdn, and Temple 
and Yarwood. As has been previously pointed out, the 
concept of a single velocity correction foi-mida is feasible in 
only two cases, namely, small stream Mach niunbers and 
vanishingly small disturbances to the main stream. It is 
desirable then to seek a single velocity correction formula 
that combines the featiu-es of these two cases. From this 
point of view, equation (36) or equation (36) is not the best 
choice, A better choice of a single velocity correction for- 



mula appears to be the following combination of equations 
(35) and (36), based on the arithmetic mean otfir) and gir): 



(37) 



In a later section, still another combination referred to as 
"the geometric-mean type of approximation" is introduced; 
in the section dealing with the Glauert-Prandtl approxima- 
tion, certain features of the foregoing arithmetic-mean type 
of approximation and of the geometric-mean type are 
discussed. 

At this point it is desirable to discuss the practical applica- 
tion of equation (37). According to equation (16), 



and 



in 



(38) 



Equation (37) then yields, for a given set of values of the 
stream Mach number Mi and the local Mach mmiber M, a 
value for the ratio (g/ffO* of the local velocity g and the 
stream velocity gi ^ incompressible fluid. Table 2 shows 
corresponding values of (g/gOc and (gJgOi for various values 
of the stream Mach mmiber Mi with 7=1.4 03=2.5). This 
tabulation is performed, for the purpose of comparison, for 
the three cases represented by equations (35), (36), and (37). 

Values of {gjQi)u {qJQi)c, and ^^^^\' > obtained from equations 

(37) and (38), are plotted against the local Mach number M 
in figure 2 for various values of the stream Mach munber Mi. 
Table 2 also shows values of the pressure coefficients Cp.o 
and Cpjfj calculated by equations (18a) and (18b) for 
these corresponding values of (gJqOt and (s/qOe- Figure 3 
shows the curves of pressure coefficients corresponding to the 
curves of velocities of figure 2: Useful cross plots of the 
curves in figure 3 are shown in figure 4, in which Op,Ui is 
plotted against Mi for various values of Cp,o. In addition, 
curves are shown in figure 4 for {Qj,,it^, and ((7p^j)„ 
calculated by equations (18c) and (18d), respectively. The 
ciurve for {Cp.u^, corresponds to the sonic value M=l or 

T=T,=^ and in effect divides the region of flow into a sub- 
sonic and a supersonic part. The curve of {OfM^m corre- 
sponds to the maximum value M= <» or r=l and represents 
the outer limit of the supersonic region (or ■& perfect vacuum) . 
In order to exhibit the main differences between the various 
correction formulas (35), (36), and (37), the ratios of the 
sonic values (C^^j), and the corresponding incompressible 
values Gfja are plotted against the stream Mach number 
Ml in figure 6. 




(b) (i2/in)^iatIootIoodTeIi>dt7tostieamTaIodt7>lnoa(mpiessIbIa 
Vtmms 3.— Oantbiiied. 
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(b) Cjfcams. 
FiGUBE S.— Oonoliided. 
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FiauBB 4.-Pressnre ODemdmt C,, j/, against stream Mach nnmber Mi with conespondine values of C^j; and linos oT constant laeal Madi nnmber Id. 



Observe in figure 2 that the (g/gOf-curres have maxi- 
mum points. This fact means that the value of {sJq_i)c 
associated with a value of (g/gi) < is not unique. Analytically, 
the criterion for the mamnum point is equivalent to 



dr 

or, from velocity correction formula (37), 

(l-r)^i-(2i3+l)T+l=0 



(39) 



For /3=2.6 this equation has only one positive root, 

or M» 1.16. It is iuteresting to note that velocity correction 
formula (36) yields as the criterion for the maxiinum point 

l-(2j3+l)r=0 
The root of this equation is T=T,=:r^^ and, for /3=2.5, is 

T=g or M=l. Velocity correction formula (35) yields no 
•maximum value of t or M. 




.a 



.6 



I.O 



FiousE 5.— The tatio of ChMi), to 0,a against li\ for the varioos apptoxlmatloiu. 



ON THE FLOW OF A COMPHESSIBLE FLtttD BY THE HODOGRAPH METHOD — 



297 




0 .1 .2 .3 .4 .5 .6 .7 .8 .3 1.0 /./ /-? 

M 

FiauBE 0.— Tbo ftmctlon V against Maoh nnmber Af for several valnes of the adlabatio Index i. 



Meaning can be given to the value T=g (M=l) in the 

case of equation (34) with reference to the original inter- 
pretation of the flow pattern as that of a somce. It can be 

shown that the acceleration ^2^^ along a streamline ^is 

infinite at all points for which the local Mach nxunber is 

unity ^'■^^ 8, flow discontinuity exists there. 

In the case of the vortex flow pattern (equation (33)), no 
flow discontinuity occurs for M<oo. The velocity correc- 
tion formula (37) suggests a "limiting" value ikf«<1.16 for a 
spu'ol flow, since equation (39) is analogous to a condition 
of infinite acceleration. Thus, the existence of a mixed 
subsonic and supersonic region of flow without discontinuities 
is indicated. Since the occurrence of this limiting value of 
Af is a consequence of the simple form assumed for the 
velocity correction formula, no imdue significance should 
be attached to any particular value at the presemt time. 

THE CHAPLYGIN APPBOmfATION 

From the point of view of the present paper, Ohaplygin's 
approximation for subsonic speeds assumes a simple and 
lucid form. Chaplygin introduces in place of g a new inde- 



pendent speed variable ij equivalent to the quantity given 
on the right-hand side of equation (33), namely, 

The hodograph flow equations (7) then assume the form 



where 



"bi, -^W-bO 



(40) 



jPf v. l-(2ff-|-l)T 

=l-i3(2|8+l)r*-|i8(2/S+l) (2^+2)7^- , . . 

Values of the function F{t), for several values of y (or /8), 
are given in table 3 and are plotted in figure 6 against the 
local Mach number M. Chaplygin noted that, ia the case 
of air (j8=2.5), F(t) differs but little from unity over about 

one-half the subsonic range O^t^ ~ His approximation in 

the range of low subsonic speeds consists in neglecting powers 
of T higher than the first or in replacing F(t) by unity. 
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Equations (40) can ihea be •written in the Cauchy-Riemann 
form 

&0~d log rj 
b log ij do 

and 4>+if therefore is an analytic function of the complex 
variable 6+i log ij. Chaplygin's approximation thus leads 
to the velocity correction formula 



1-4^1 



(41) 



wherepowers of r higher than the first are neglected through- 
out. The use of equation (34) instead of equation (33) also 
leads to this result to the same order of approximation. 

THE TON KIBMIn APFROXIMATION 

Von K&rmdn's approximation corresponds to the case 

Y=— 1 (oT P=—^- It follows at once from the integral 

expressions for /(t) and gir) given by equations (26) and 
(27), respectively, that for this case 

or, with the use of equation (16), 

Ar)=gir) log i 

This function, plotted against M, is included in figure 1(b). 
Corresponding to equations (35) and (36), there is a single 
equation 

(l)r(l)ci+(i-ry'' 

Replacing r by n and n by according to equa- 

tion (16) yields 

1+(1-Mi')'« 



=(^) 



(42) 



Then, by solving for (g/gOe in terms of (j/gi)* and the stream 
Mach number Mi, 



where 



^-nT(i=M?)^ 



(43) 



The pressure coefficient Cp^^, expressed in terms of the 
incompressible pressure coefficient Cp,o, is easily obtained 



from the general formula (18b) by putting 7= — ! and 
making use of equations (43) and (18a). Thus, 



(44) 



Observe that for this case the fimction F{t) introduced by 
Chaplygin and given in equation (40) is exactly equal to 
unity. From the point of view of the present paper then, 
von Eton&n's approximation appears to be equivalent to 
that of Chaplygin, who approximates F(t) by unity. It 
follows that the range of validity of von Ktonfin's approxi- 
mation and that of Chaplygin, in a strict sense, coincide, 
furthermore, it is pointed out that the von K&nn&n approxi- 
mation does not permit a supersonic region. Von K&rmfi,n's 
choice of y= — 1 has the advantage, however, of yielding 
simple explicit expressions for (g/gi)« in terms of (qJgOt and 
for C7pjifi in terms of <7p,o. Several values of Of,ui calcu- 
lated by equation (44) are included in figure 4. For the 
purpose of comparison with the other approximations, there 
is plotted in figure 5 the ratio of (Cpjui), to (7p,o against the 
stream Macb number Mi in the case of von Ktan&n's approxi- 
mation. The values of (7p,o are obtained with the use of 
velocity correction formula (42) for the local Mach number 
M=l, but the values of (Cp,ui)$ are calculated AVith 7=1.4. 

THE TEMPLE-YABWOOD APPROXIMATION 

The functions ^ and ^ related by the first-order simultane- 
ous equations (21) separately satisfy the second-order 
equations 



(46) 



In terms of the nondimensional speed variable r and with 
the values of \i(g) and Xa(2) for the adiabatic case given by 
equations (25), tihese equations take the form 



1 (l-ry i^<}> b 
4 T be*"'" 5 



) r t(I-t)'^^ b<f] 
;Ll-(2/3+l)r brj^ 

1 1-(2/3-H)t . a r r bf]_^ 
4 t(1-t)<«-i dO*'^bTL(l-T)'' &rj~"/ 



(46) 



Formal solutions of these equations were given by Chaplygin 
in the form of two infinite series 



^=5<?+SS„^„(r) sin (m0+t„) 



<^=-5^(t)-SB»*«(t) cos ime+e„) 



(47) 



where the functions ^mfr) and <l>mir) are obtained from 
hypergeometric series and B, B„, and are arbitrary 
constants. 

A disadvantage of the formal solution, as remarked by 
Temple and Yarwood, is that it is imsuitable for numerical 
computation because the hypergeometric functions involved 
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are complicated and are not tabulated. Temple and Tar- 
wood therefore looked for approximations that are of 
practical value in calculations of compressible flows. By 
means of a skillful analysis, they found such approximations 
and showed that the simplest forms for ^„ and <l>„ are of 
the type 

^«(r)«h(r)l» 

(48) 

^(t) «logf(r), 
where ij(t) and {(t), independent of the index m, are 

,=f=(l-|r)2 (49) 

Significantly, from the point of view of the analysis of the 
present paper, the functions ij and ^, approximated by 

^1— 2 are none other than the functions defined on the 

right-hand sides of equations (33) and (34). The approxima- 
tion of Temple and Tarwood then leads to the same velocity 
correction relation as was obtained by means of Chaplygin's 
approximation (equation (41)). 

The velocity and pressure-coefficient correction formulas 
obtained by Temple and Tarwood are more involved than 
the explicit expressions (43) and (44) obtained by von 

K&em&a. Replacing t in equation (41) by '"i^^^^ thus yields 



X—-r Ti 



(50) 



where 



_ Ml' 



The solution of this cubic equation for (s/gi)e is 



COSg (tT+O-) 



cos o- 



(51) 



where 



COS 0-= 



(-f-)(l")'"(D. 



and 0<<r^^' The pressure coefficient Gj,^^ is then calcu- 
lated by equation (18b). Some values of the pressure 
coefficient O^u, calculated with the aid of equation (51) are 

(O \ 

shown in figure 4: a curve of /y "' plotted against Mi is 

included in figure 5. It is remarked that, with the use of 
equation (39), the velocity correction formula (50) yields a 
limiting value M*» 1.35. 

APPROXIMATION BASED ON GEOMETBIC MEAN OF dL AND d2 

Without going into its deep significance in the present 
paper, it is of interest to introduce another function related 



to L and L and to the general particular solutions. This 
fimction, which like L and L reduces to log j for t=0, is 
defined by 

H(T)=J(dLdLy'^ (52) 

It is remarked that S(t) is closely related to a function K{t) 
employed by Temple and Tarwood (reference 7) in the 
determination of their approximation. In the next section, 
it will be seen that the function H(t) plays an important role 
in connection with the Prandtl-GIauert approximation. 
From equations (26) and (27), 



and 



Then, 



and, from equation (52), 

F(T)=log2+A(T) 



(63) 



(54) 



where 



The function A(r) can be obtained in a closed form for any 
value of 7 (or jS) and is 



AC 



r(l-r)»^-f-(l-jy'n[(l-r)^/*-(r.-r)''»]^ 

y) — log L ^ i 



where r,= 



2/5+1 



(55a) 



and where this expression is valid in the 

M* 



subsonic range O^t^t,. With r replaced by 20+]ji^^^ 
O^M^ 1, the expression for A(t) becomes 



r) — log — 2 -dz 1-v^: 



i+V^i,-i+V^(i-W^ 



log 



(56b) 



It is observed that, for the supersonic region t, ^rgl or 
M>1, H(t) as defined by equation (52) becomes a complex 
function; but, for present purposes, only the real function 
of the subsonic range is utilized. 

The function H(r) may be utilized to obtain a velocity 
correction formula in the same manner as the functions 
Z(t) and L{t). Thus, analogous to equation (35), (36), 
or (37), 
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It is instructive to compare equation (56) with the approxi- 
mation given by equation (37). Equation (37) may be 
written as 

and equation (56) may be written as 



(l)r 



Thus, the power of the exponential is in one case the integral 
of the arithmetic mean ^^^"^ and in the other case the in- 
tegral of the geometric mean (dL dL)^^. Table 1 shows val- 
of the fimctions -^^^^^i^^ and A(t) in the case of 



ues 



au: 



^7=1.4, /3=2.5, and t.=^^ and figures 1(a) and 1(b) show 

these functions plotted againstr andM, respectively. Observe 
that these fimctions, and consequently the velocity correction 
formulas (37) and (56), differ only slightly-in the subsonic 
range 0<M<1. Figure 5 exhibits graphically a comparison 
of the velocity correction formulas (37) and (56) for M=l. 
The limiting value of M (defined by equation (39)) is M=l 
in the case of equation (56) as compared with M«s 1.16 in the 
case of equation (37). 

COMPABISON OF RESULTS OF PRESENT PAPER WITH PRANDTL- 

GLAUERT APPROXIMATION 

The well-'known PrandtiL-Glauert approximation is based 
on the assumption of vanishingly small disturbances to the 
main stream. The Prandtl-Glauert velocity correction 
formula may be expressed as 



\ gl Jc _ 1 



(57) 



where g— gi is vanishingly small. The left-hand side of this 

equation is actually the differential coefficient 

evaluated at the main stream velodty g=2i (or t=ti). An 
exact form of the Prandtl-Glauert approximation then is 



R(g/gl)c 1 ■ 1 



(58) 



The differential coefficient in equation (58) is now evaluated 
for the various approximations treated in the present paper. 

For the arithmetic-mean approximation of the present 
paper given by equation (37) (7 or |3 arbitrary), 

r ^(g/gi)c 1 2 

U(2/gi)J-r. ■ 1-(2/3-1-1)ti 



i+(i-m4i+^) 



=i+|m*+|m^-I-^m,« 



^128 j3* ^ ^ 



(59) 



For the Chaplygin or the Temple-Yarwood approximation 
given by equation (41) (7=1.4 or jS=2.5), 



[ dm). -] _ 

U(g/gl)J r^i 



15 



11 TI^S 



(60) 



For the von Kfirmfin approximation given by equation (42) 
(7=-l or p i) 

r d(qjqi)c ~\ ^/i__-vi/2 



(61) 



For the geometric-mean approximation of the present 
paper given by equation (56) (7 or /5 arbitrary), 



U(g/gl)Jr-r, L1-(2/3+1)tJ 



1/2 



(62) 



Equation (62) is independent of the value of the adiabatic 
iudex 7 and includes the von KArmdn approximation. Ob- 
serve that the geometric-mean approximation yields the 
Prandtl-Glauert result exactly, whereas the arithmetic-mean 
approximation yields the Prandtl-Glauert result insofar as 
terms inclusive of Mi' are concerned. The Chaplygin or the 
Temple-Yarwood approximation contains the Prandtl- 
Glauert result only insofar as the Afi'-term is concerned. 

KfiSUMfilAND CONCLUDING REMARKS 

1. Basic elementaay solutions of the hodograph equations 
have been employed to provide a basis for comparison, in 
the form of velocity correction form\ilas, of corresponding 
compressible and incompressible flows. 

2. The velocity correction formulas obtained by Ohaplygm, 
by von K^rmdn, and by Temple and Yarwood have been 
unified by means of these basic solutions and shown to be 
essentially equivalent. 
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3. In the present paper two types of approximations have 
been introduced by means of the basic elementary solutions, 
namely, the "arithmetic-mean" type and the "geometric- 
mean" type. These approximations include those obtained 
by Chaplygin, by von Ktondn, and by Temple and Tarwood. 

4. The approximations discussed in the present paper 
have been compared with the well-known results of Prandtl 
and Glauert. For this purpose, it has been emphasized that 
the Prandtl-Glauert result is valid for vauishingly smaU dis- 
turbances and, in a strict sense, is the slope term in a Taylor 
expansion in a quantity which measures the disturbance. 
It was found that the arithmetic-mean type yields the 
Prandtl-Glauert result to a higher order of approximation 
than the Ohaplygin or the Temple-Yaiwood type and that 
the geometric-mean type contains the Prandtl-Glauert result 
exactly. The two types of approximations introduced in 
the present paper then appear to be preferable to the others 
as a basis for extrapolation into the range of high stream 
Mach numbers and large disturbances to the main stream. 

6. The results of the present paper have been obtained 
mthout consideration of any particular boundary. The 
actual boundary problem of determining the flow past a 
prescribed body is of a high order of difficulty and involves 
in general all the particular solutions of the hodograph 
equations. 

6. The particular solutions discussed in the present paper 
are well-behaved functions in both the subsonic and the 
supersonic regions. The hodograph equations give no reason, 
in general, to suppose that a discontinuity necessarily occurs 
in the solution when local sound speed is attained. Rather, 
it appears that the first breakdown of the solution is as- 
sociated with the vanishing of the Jacobion of the trans- 
formation from the physical to the hodograph variables. 
Indeed, von Kdrmdn has made an equivalent suggestion ia 
that the appearance of infinite accelerations in the flow 
solution is a condition for flow discontinuities. Interesting 
speculations on this matter are suggested by the results of the 
present paper since the "limiting" curves discussed in the 
present paper are defined by a condition that is equivalent 
to the condition for infinite acceleration. The arithmetic- 
mean type of approximation thus yields a limiting value of 
the local Madi number iWasl.16, and the geometric-mean 
type of approximation yields a limiting value of the local 
Mach number M=l. The value M=l appears to be exact 
for vanishingly small disturbances; that is, local Mach 
number M=8tream Mach number Mi=l (Prandti-Glauert 
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approximation). However, for finite disturbances to the 
main flow due to the presence of a body in the fluid, infinite 
accelerations may occm-, for stream Mach numbecre less than 
imity, in regions where the local Mach number is greater 
than unity. In this regard, the arithmetic-mean type of 
approximation, considered as an extension of the Prandti- 
Glauert relation to finite distiu-bances, indicates the pos- 
sibility of a mixed subsonic and supersonic flow without 
discontinuities. It is important, however, to recognize 
that in general the limiting value of the local Mach nmnber 
M is a function of shape parameters and is a result of the 
blending of many particular solutions of the hodograph flow 
equations according to the boundary conditions. 
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TABLE 1 



.—VALUES OF /, ff,^. h, AND THEIR EXPONENTIALS FOR 7=1.4 



M 



.1 

.3 

.3 

.* 

.S 

.6 

.65 

.TO 

.75 

.SO 
.825 
.850 
.876 
.900 
.925 
.950 
.96 
.98 
LOO 

1.02 
1.M 
LOS 
1.03 
LIO 
1.12 
LIS 
LIS 
LIO 
1.30 

1.40 
LM 
3.00 
2.50 
3.00 
4.00 
ILOO 



.00200 
.007M 
.01768 
.03101 
.0<782 
.06716 
.07792 
.08925 
.10112 

.U3*8 
.U982 
.12626 
.13279 
.13942 
.14612 
.15290 
.15563 
.16113 
.16667 

.17224 
.17785 
.18349 
.18915 
.19485 
.20056 
.20917 
.21782 



.31034 
155566 



.761900 
.8331 
LOO 



0 

-.00250 



.03196 
.03831 
-.06847 
.08186 
.09457 
-.10787 
-.12167 

-.13688 
-. 14313 
-.15046 
-.16785 
-.16530 
-.17281 
-.18038 



.18946 
.19556 

.20168 
.20778 
.21391 
.22003 
.22616 



-.34144 
.25059 
.25665 
-.28673 

.31604 
.34438 
.46775 



-.62170 



-.74934 
-.84019 



0 

-.00251 
—.01013 
-.02316 
-.04208 
-.06760 
-.10059 
—.12023 
-. 14216 
-.16667 

-.19363 
-.20822 
-.22354 
-.23964 
-.25652 
-.27423 
-.29280 
-.30047 
-.31624 



-.34958 
-.36718 
-.38642 
-.40432 
-.42390 
-.44117 
-.47504 



—.81292 
-.99030 
-2.39742 
-5l 12014 
-9.99177 
—3L 27238 
-8a 81740 



2 



0 

-.00251 
-.01001 
-.02256 
-.04020 
—.06301 
— . 09123 
-.10740 
-.12502 
-. 14412 

-.16476 
-.17568 
-.18700 
-.19875 
-.21091 
—.22363 



-.34193 



-.26409 

-.27662 
-.28748 
-.29967 
—.31218 



—.37999 
-.39464 
-.47406 

—.86448 
—.66734 
-L«259 
-2.83970 
-6.30834 



—40.78337 



0 

—.00250 
-.01001 
-.02266 
-.04020 
-.06306 
-.09133 
-.10788 
-.12511 
-.14484 

-.16605 
-.17740 
-.18927 
-.20173 
-.21487 
—.22876 



-.21976 



—.27787 







e > 


LOO 


LOO 


1.00 


.99750 


.99749 


.99750 


. &9016 




99001 


.97828 


! 97711 


;9T768 


.96241 


.95879 


.96060 


.94321 


.93463 


.93890 


.92140 


.90430 


.91281 


.90977 


.88671 


.89817 


.89775 


.86748 


.8SMS 


.88544 


.84666 


.88578 


.87296 


.82397 


.84810 


.88664 


.81202 


.83889 


.86032 


.79968 


.82944 


.85398 


.78681 


,81976 


.84764 


.77374 


.80985 


.84130 


.76016 


.79969 


.83497 


.74617 


.78933 


.83315 


.74047 


.78511 


.82741 


.72888 


.77658 


.82237 


.71705 


.76791 


.81737 


.70498 


.76910 


.81239 


.69268 


.76015 


.80742 


.68017 


.74106 


.80250 


.66743 


.73186 


.79768 


.66449 


.72251 


.79273 


.64136 


.71303 


.78660 


.62130 


.69869 


.77834 


.60087 


.68387 


.77364 


.68700 


.67392 


.76072 


.61613 


.62»7 


.72903 


.44366 


.66886 


.70866 


.3n47 


.61307 


.62641 


.09096 


.23869 


.6n64 


.00698 


.06844 


.53642 


.00004 


.00496 


.49393 


.00000 


.00000 


.47268 


.00000 


.00000 


.43162 


0 


0 
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TABLE 2.— TALUE8 OF (a/gi)., (g/ji),, -^j-', Cp.„, AND C,.j„ 
FOR 7=1.4 AND FOR VARIOUS VALUES OF M, 



M 


a2 


as 


0.* 


as 


a££ 


as 


a7 


a8 


a9 


LO 


LI 


L2 


T 


aooTU 


0.017SS 


Oi 03101 


a 04782 


aosTOs 


aosTie 


a08926 


ail34S 


a 13942 


a 16687 


a 10135 


a22360 










LOO 


L 49262 


LS7eeo 


2.44S48 


2.68106 


2.00907 


3.36349 


3.78124 


4.19121 


4.68266 


4.06184 


6.30790 




Eq.(3Q 


1.00 


L 47471 


LS3120 


2.33332 


2L 62518 


Z7D706 


3.04061 


3.3S365 


3.6ST93 


3.S0599 


3.99122 


4.14722 




Eq.(38) 


LOO 


L 47331 


L 91444 


2.31288 


2.49242 


2.65746 


2.93870 


3.14736 


3.27693 


a 31937 


3.27690 


3.14780 




Eq.(37) 


LOO 


L 47401 


L 91782 


Z32296 


Z 60376 


2.68213 


2.98917 


3.23912 


3.42841 


3.66441 


3.61691 


3.61300 






LOO 


L 01263 


L03066 


L06446 


L0686S 


L 08461 


L 12188 


L 16737 


L 22210 


L2S926 


L37D2S 


L46907 




Eq.(37) 


0 


-L 17271 


-Z 67803 


-4.39614 


-6.29383 


-a 10382 


-7.93614 


-9.49190 


-la 76400 


-IL 63383 


-IZ 07602 


-IZ 05442 


(Eq, (18b)) 




0 


-L 21260 


-2.82393 


-4.76600 


-a 81393 


-S. 92214 


-a238S3 


-IL 62214 


-14.000n 


-16.31357 


-18.51393 


-2a 67000 




(«/?! 




a6a997 


LOO 


L 34426 


1.04107 


L 79621 


L948S8 


2.24672 


2.63331 


Z80797 


3.07016 


3.31968 


3.66612 




Eq.(36) 


a 67810 


LOO 


L30277 


L6S224 


L 71233 


1.83fi67 


2.06177 




Z4329g 


Z6S0S7 


170646 


Z81226 


(J/fl)( 


Eq. (36) 


a 67876 


LOO 


L29942 


L 66971 


L 69172 


1.80374 


L99463 


Z 13628 


Z22366 


Z26304 


Z22362 


Z 13866 




Eq. (37) 


0.67843 


LOO 


L 30110 


L 67696 


L70200 


L 81966 


Z02793 


Z 19762 


Z32693 


Z4U40 


Z 46317 


Z 46123 


(?/?l). 




a837£3 


LOO 


L 01779 


L 04131 


L065S6 


L 07107 


L 10789 


L 16280 


L20725 


L 27319 


L 35318 


L 46076 


(Eqfjito)) 


Eq. (37) 


0.S3S73 


0 


-a 69286 


-L 48365 


-L 89880 


-2.81113 


-3.11260 


-3.82009 


-4.40995 


-4.814S5 


-5.01804 


-6.00863 


(Bq. (18b)) 




0.6£7si 


0 


-a 74111 


-LS2984 


-2.11883 


-Z 62651 


-8.69238 


—4.78821 


-6.88302 


-6.94746 


-7.96968 


-8.90571 


Afi-0.4 






a 76614 


1.00 


L 23934 


L36S40 


L 47174 


L69669 


L 91800 


Z12040 


Z 31840 


Z 50674 


Z68536 




Eq. (36) 


0.52051 


a 76769 


1.00 


1. 21461 


L 81438 


L40903 


L58260 


L 73510 


L 86764 


L 98105 


Z 07744 


Z 16865 


(j/»l)< 


Eq.(3fl) 


ae2235 


a76967 


LOO 


L 20801 


L 30190 


LS88I1 


1.63603 


L 64401 


L 71116 


L73387 


L71U6 


L04422 




Eq.(37) 


0.S3143 


0.76S67 


LOO 


L 21124 


1.30813 


L39863 


L 66863 


L68897 


L78766 


1.S633S 


L8S512 


L8SS91 


(?/ffi). 
(«/ji)( 




a97025 


0.98263 


LOO 


L023LJ 


L03690 


L09236 


L 08862 


L 13261 


L 18614 


L250g2 


L32961 


L42540 


(Eq. (18a)) 


Eq. (37) 


a 72811 


a40930 


0 


-a 46710 


-a 71120 


-a 96589 


-L 42030 


-L 86262 


-Z 19669 


-Z 43491 


-Z 66181 


-Z 51923 


(Eq. (18b)) 




a 76643 


a 43714 


0 


-aE3420 


-a 81188 


-L 11269 


-L 74162 


-Z 38866 


-3.03429 


-3.66205 


-4.2593S 


-4.81750 




(?/jO. 


a4082« 


a60936 


a80696 


LOO 


L09461 


L18763 


L3690e 


LM370 


L7U06 


L87034 


Z02281 


Z 16691 




Eq.(3«) 


0,428*7 


a6S202 


a82338 


LOO 


L08222 


L 16017 


L30307 


L42871 


L6376S 


Le3U5 


L 71051 


L77737 


(f/!l)f 


• Eq.(38) 


a43340 


0.63706 


a82781 


LOO 


L07773 


L 14909 


L27069 


L36092 


L 41662 


L 43631 


L 41661 


L 36110 




Eq.(37) 


a 43049 


a634(4 


aS2660 


LOO 


La790S 


L 16462 


L28679 


L39440 


L 47688 


L 63013 


L66662 


L66539 


(«/9i). 
(«/li)< 




0.04834 


0.96032 


a 97741 


LOO 


L 01348 


L02S6g 


L0e393 


L 10707 


L 16936 


L22267 


L29941 


L 39318 


(Eq. (I8a)) 


Eq.(37) 


0.81468 


a69736 


a 31838 


0 


-0.16636 


-a 33316 


—0.66683 


-a 04436 


-L 17822 


-L 34127 


-L 42307 


— L-41924 


(Eq. 08b)) 




a 87771 


a66366 


a36S46 


0 


-0.19664 


-a 40000 


-a 82766 


-L 26764 


-L 70667 


-Z 13343 


-Z 63960 


-Z 91903 


JW1-O.6S 


(1/30 


• 


a37299 


a66672 


a7S726 


a 91363 


LOO 


L 08604 


L 26081 


L 41036 


L6e328 


L70923 


L 84810 , 


L97978 




Eq. (35) 


a 39601 


a68399 


a 76083 


a92404 


LOO 


L07202 


L20408 


L 82017 


L42aS8 


L6D722 


L5S067 


Le4235 


{g/«i)< 


Eq.(38) 


a 40122 


a 69110 


a 76811 


a927S8 


1.00 


L06621 


L 17906 


L 26278 


L 31434 


L 33178 


L 31436 


L262S1 




Eq. (37) 


0.39861 


a 68763 


a 76446 


a92696 


LOO 


L 06011 


L 19160 


L 29114 


L36S67 


L 41680 


L 44133 


L 44010 


mot 




a93573 


a 94766 


0.96441 


a9S669 


LOO 


L 01490 


L 04978 


L0I234 


L 14393 


L2064a 


1,28222 


1.37467 


(Eq. (18a)) 


Eq.(37) 


a841U 


0.66481 


a 41660 


a 14283 


0 


-a 14300 


-a 41967 


-0.66704 


-0.88751 


—1.00732 


-L 07743 


-L 07418 


(Eq. (18b)) 




a 91836 


a72686 


a 47^9 


a 16737 


0 


-a 17497 


-0.64078 


-a 91731 


-L 29299 


-L 66828 


-Z 00690 


-Z33068 
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TABLE 2.— Continued 



M 


a4 


a« 


ass 


a7 


ttTS 


as 


ass 


ao 


LO 


1.1 


L2 


L3 


T 


a 03101 


a 06716 


ao77a3 


aOS926 


a 10112 


ausis 


am26 


a 13942 


416667 


419485 


422380 


425203 




te/5i). 


0.67917 


LOO 


L07707 


L 16277 


L227D1 


L 29981 


L 37107 


L44093 


1.57527 


1.70324 


1.82400 


1.93939 




Eq.(30 


0L70B71 


LOO 


L0634S 


L 12318 


L 17914 


L 23146 


L 28018 


L32669 


L40590 


1.47438 


1.63201 


1.68014 




Eq.(38) 


0.72011 


LOO 


L 05612 


L 10633 


L 14860 


LIS^ 


L 21245 


L23289 


1.24008 


L23272 


1.18461 


1. 10091 




Eq.(37) 


0. 71601 


LOO 


L0S979 


L 11447 


L 16382 


L 20767 


L24S85 


L27841 


1.32521 


L 34816 


1.34703 


1.32363 






0.9S02S 


LOO 


L 01631 


La3437 


L05432 


L 07630 


L 10051 


1.12713 


L18869 


L26339 


L36449 


1.40042 


c,» 

(Eq. (18a)) 


Eq.(37) 


0.1SS72 


0 


—a 12315 


—0. 24204 


-a 36448 


-a 45847 


—0.55214 


—a 63433 


—4 76618 


—4 81761 


—481402 


-0. 7490B 


(Eq. (18b)) 




aseiSQ 


0 


-0.15786 


-a 31929 


-a 48310 


-a 64771 


-a 81233 


-497599 


-L 29437 


-1.69763 


-L 88099 


-2. 14161 














M°0.65 












• 


(S/«i). 


aosoM 


a92S43 


LOO 


L07029 


L 13924 


L206S1 


L 27296 


L33766 


L4S266 


L 68138 


L 69405 


L 80000 




Eq.(35) 




a04O3O 


LOO 


L0561B 


L 10878 


L 16797 


L20377 


L 24630 


L32201 


L33638 


L 44067 


L 48583 


(2/}i)i 


Eq.(38) 


a 68212 


a946S5 


LOO 


LO4708 


L08786 


L 12143 


L 14802 


L 16723 


1.18271 


L 16723 


L 12168 


L 04808 




Eq.pT) 


a 67169 


a04356 


LOO 


L 05159 


L09S16 


L 13953 


1.17555 


L 20612 


L25044 


L27209 


L 27110 


L 24789 


(j/jiJT 




a 03501 


a9S396 


LOO 


L 01778 


L 03741 


L05904 


L0S288 


L 10909 


L 16963 


L 24314 


L 33274 


L442S3 


(Eq. (18a)) 


Eq.(37) 


aM179 


a 10969 


0 


-a 10584 


-a 20996 


-a 29853 


-a 38192 


-a 45473 


-466380 


-461821 


-461670 


-466723 


{Bq. {18b)) 




a 64071 


a 14003 


Q 


-a 14333 


-a 28882 


-0.43476 


-assose 


-a72S71 


-L 00862 


-L 27767 


-L 62913 


-L 76024 


Afi=0.7 


ft/fi). 


aSS912 


a 86747 


ae3433 


LOO 


L06442 


L 12756 


LIS936 


L 24980 


L 38650 


L 47763 


L 68280 


L0S235 




Eq.(35) 


a 63186 


aS9032 


a946S3 


LOO 


L 01983 


L09641 


L 13078 


L180a3 


1.25176 


L 31208 


L 30398 


1.40081 


(S/Ji)i 


£q.{38) 


a 65146 


agoi29 


a05604 


LOO 


L 03876 


L07099 


L09640 


1. 11476 


L12S63 


L 11476 


L 07114 


L 00090 




Eq. (37) 


aeiuo 


a8972S 


a95094 


LOO 


L04429 


L 06362 


L 11788 


L 14696 


L 18900 


L20960 


L 20374 


1. 18600 


(«/?i). 
(ffAi)< 




a91S69 


a966T8 


a9S253 


LOO 


L 01928 


LO4055 


L 06394 


L08967 


L 14920 


L 22140 


L 30940 


L 41772 


(Eq. (18a)) 


Eq.(37) 


0158836 


a 19489 


ao95n 


0 


-a 00054 


-0.17423 


-a 24966 


-431649 


-411394 


-440336 


-446105 


-440810 


(Eq. aSb)) 




aTosso 


a 35516 


aisooi 


0 


-a 13082 


-0.26264 


-a 39397 


-a 62140 


-477907 


-402140 


-L 24781 


-L 46680 


Jvri-o.7s 


(»/Jl)« 


a 55374 


a 81497 


0.87778 


a 93947 


LOO 


L 05931 


L 11737 


L 17116 


L2S380 


L 38810 


L 48701 


L58054 




Eq.OS) 


a601S7 


aS4S07 


a 90101 


a95253 


LOO 


L04438 


L08668 


L12404 


L 19235 


L250?f> 


L29026 


1.34000 


(!/«i)« 


Eq.(38) 


aa27i6 


asTosa 


a 91942 


a9626S 


LOO 


L 03105 


L06650 


L 07317 


LOS740 


L 07317 


L 03120 


490362 




Eq.(37) 


a 61439 


a 85924 


a 91062 


a 95759 


LOO 


L 03768 


L 07047 


L 09831 


L13888 


L 15840 


L 16749 


1.13030 


(!7/fl). 
(?/fi)f 




aoom 


aoisis 


a96394 


aosios 


LOO 


L020S4 


LW381 


L08906 


L 12745 


L1C329 


L2840S 


L390SS 


(Eq. (IBa)) 


Eq.(37) 


a62252 


a 28171 


a 17077 


aos302 


0 


-0.07678 


-a 14691 


-a 20628 


-429659 


-434189 


-433978 


-420131 


(Eq. (18b)) 




a76361 


a 35197 


a 23700 


a 11937 


0 


-a 12005 


-a 23997 


-435893 


-459124 


-481227 


-L 01872 


-L 20863 


Ari-o.s 


(!/3i) 


« 


aS2Z74 


a76S34 


aS2S84 


a8S6S7 


a94402 


LOO 


L 05482 


. L10642 


L 21193 


L 31038 


1.40376 


L 49201 




Eq.(35) 


a«7B3I 


aS1205 


aSS359 


a 91207 


aS57«3 


LOO 


L03a56 


1.07629 


L 14171 


L 19725 


L 24400 


L 28312 


(9/!i)( 


Eq.(36) 


a60327 


a84435 


a 89172 


a 03371 


a96990 


1.00 


L023i2 


L01036 


L06467 


LO1035 


L 00014 


4 93400 




Eq.(37) 


0lS9209 


aS2806 


a 87766 


a022S3 


a 06370 


LOO 


L 03162 


L05S43 


L09734 


L 11633 


1.11645 


L 09510 


(j/?i). 




a88287 


a 92910 


0.94425 


a 96103 


a97958 


LOO 


1.02249 


L 04723 


L 10443 


L 17383 


L2SS40 


1.30247 


(Bq.{18a)) 


JBq.(37) 


a64943 


a 31433 


a229S9 


ai483S 


a 07128 


0 


-a 06424 


-412027 


-a2M16 


-424619 


-a 24423 


-419024 


C,jf, 

(Eq. {18b)) 




a81S20 


a43547 


a32940 


a22083 


a 11071 


0 


-a 11065 


-a 22036 


-a 43464 


-a 63867 


-482904 


-L 00411 
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TABLE 2.— Concluded 



M 


a4 


16 


1825 


186 


1876 


19 


1925 


196 


LO 


LI 


L2 


L3 


T 


a 03101 


aoene 


1II982 


112628 


113279 


113942 


114612 


116663 


116667 


110485 


122360 


125262 






a«)872 


174871 


LOO 


L02663 


1.05276 


L078n) 


L 10134 


L13972 


L 17942 


L 27624 


L 36610 


L 46203 




Eq. (36) 


a664B4 


0.70602 


LOO 


L 01006 


L 03738 


L06606 


L07205 


L 09476 


L 11917 


L 17363 


L 21950 


L25781 




Eq.(36) 


a60068 


183380 


LOO 


L 01093 


L 02021 


L 02786 


L03380 


L03929 


L 04147 


L02784 


1987«t 


102292 




Eq.(37) 


a«S2fi3 


181469 


LOO 


L 01497 


L 02874 


L 04136 


L05273 


L06665 


L07962 


L09SS3 


L 09746 


L 07744 


Mil" 

(llQi)i 




a87329 


191901 


LOO 


L 01139 


L02335 


L03588 


L 04002 


L 06850 


L 09244 


L 16107 


L 21478 


L 34767 


(Eq. (18a)) 


Eq.(37) 


0.66066 


133628 


0 


-103016 


-Oi 05831 


-108443 


-110821 


-113774 


-116568 


-120633 


-120142 


-116088 


(Eq. (18b» 




o.S3m 


147330 


0 


-106329 


-110631 


-115910 


-121140 


-123409 


-138571 


-1662^ 


-174509 


-191487 






(J/Ji). 


IH9SS3 


172936 


197416 


LOO 


L02656 


L050S2 


L 07680 


L 11026 


L 14894 


L2422S 


L3308a 


L 41461 




Eq. (3£) 




178114 


198132 


LOO 


L 01799 


L03632 


1.05202 


L 07430 


L09828 


L 16171 


L 10672 


1.23432 




Eq.(3S) 


asMis 


182478 


198920 


LOO 


L 00917 


L 01673 


1.02284 


L02805 


L03023 


L 01673 


107697 


1 91205 




Eq.(37) 


aS7301 


180287 


198626 


LOO 


L 01357 


L02699 


L 03721 


L 05091 


L063n 


L08212 


L03I28 


L 06155 






a 86347 


100SS7 


198373 


LOO 


L0U82 


L0»20 


L 03721 


L 05647 


L 08013 


L 14801 


L230i<I 


L33219 


CBq. (18a)) 


Eq.(37) 


a67059 


136672 


102926 


0 


—102732 


-105266 


-107580 


-110441 


-1 13148 


-117098 


-116917 


-1I26S9 


(Eq. (18b)) 




0.86288 


ISOSOS 


106147 


0 


-105129 


-110226 


-116294 


-122305 


-130199 


-149202 


-166952 


—183274 














M 


-0.87S 














(?/«i)« 


a48322 


1 71118 


194038 


107603 


LOO 


L02463 


L048S8 


L0S269 


L 12031 


L 21133 


L 29763 


1.37926 




Eq.(3IS) 


aM468 


176734 


196396 


0.98232 


1.00 


L 01703 


L 03340 


L05630 


L07885 


L 13134 


L 17655 


L 21248 


{titiii 


Eq.(30) 


0.688T7 


181727 


198020 


199091 


1.00 


L 00748 


L 01333 


L 01870 


LQ2a86 


L 00740 


196807 


1S0164 




Eq.(37) 


a68626 


179192 


19T207 


198661 


LOO 


L 01225 


L 02331 


L03684 


L 04947 


L06764 


L06S78 


L 04733 


W/«i)< 




0.85336 


189306 


197717 


198831 


LOO 


L 01223 


L02S09 


L 01412 


L 06750 


L 13459 


L 21640 


L 31693 


(Eq. (18a)) 


Eq.(37) 


a 67936 


137288 


106608 


102660 


0 


-0.03465 


-101716 


-107601 


-110139 


-113988 


-113802 


-109690 


(Eq. (ISb)) 




a88Ml 


1S4277 


109960 


i04gn 


0 


-101943 


-109846 


-1 16642 


-0.21280 


-142701 


-169893 


-175712 


M-0.9 


(j/«i). 


147160 


160409 


192705 


196164 


197595 


LOO 


L02376 


L 05667 


L09338 


L 18220 


L28ei4 


L34609 




Eq.(3«) 


aS3646 


175449 


191783 


0.96587 


198325 


LOO 


L 01610 


L 03761 


L06079 


L1I340 


LI55SS 


L 10218 


(?/»0( 


Eq.(S6) 


0.58439 


181121 


197291 


198365 


199257 


LOO 


L00679 


LOUIS 


L 01327 


LOOOOO 


196089 


189792 




Eq.(37) 


165939 


178233 


196029 


197466 


198788 


LOO 


L 01091 


L 02129 


L0367B 


L 06170 


L053S7 


L 03463 






0.84306 


188721 


196639 


197638 


198792 


LOO 


L 01271 


L 03161 


L 05461 


L 12080 


L 20170 


1.30104 


C,.o 
(Eq. (18a)) 


Eq.(37) 


168708 


138796 


107784 


105004 


102109 


0 


-102194 


-101917 


-107487 


-111239 


-lllOM 


-107018 


(Eq. (18b)) 




0.90787 


167490 


114460 


109617 


0.04799 


0 


-1047(12 


-0.11385 


-a 18788 


-136660 


-163382 


-168711 














All 


"0.925 














(?/«!)• 


146066 


167798 


100562 


192965 


196331 


1:0.07678 


LOO 


L03205 


L06800 


L 16478 


L23704 


1.31486 




Eq.(35) 


162697 


174263 


193279 


195067 


196768 


198416 


LOO 


L02U0 


L 04397 


L 09476 


L 13766 


L 17329 


(?/«l)< 


Eq.(38) 


158103 


1806« 


196729 


0.97787 


198688 


199123 


LOO 


L00532 


L 00743 


199134 


195534 


189276 




Eq.(37) 


165336 


177387 


194991 


196418 


197721 


198919 


LOO 


L 01323 


L02556 


L 01332 


L 04219 


L 02347 


m,). 

mi)i 




183249 


187609 


195327 


19&U3 


197651 


198746 


LOO 


L 01857 


L 01138 


L 10681 


L 18662 


L284n 


C,.o 
(Eq. (18a)) 


Eq.(37) 


169380 


140113 


109767 


107046 


101506 


102160 


0 


-0.02664 


-1 05177 


-108862 


-108679 


-1 01749 


C,,i/, 

(Eq. (18b)) 




102063 


160574 


118708 


113901 


109301 


101633 


0 


-106426 


-113642 


-131038 


-147285 


-162222 
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TABLE 3.— VALUES OF F{t) FOR SEVERAL VALUES OF y 







M 


•r-1.4 
Adbbatto 


Isothennal 


r-2 

OJ-1) 
Hydraulie 
analogy 


Limiting 
incompres- 
sible! 


0 
.2 
.4 
.8 
.65 
.70 
.75 
.SO 
.85 
.90 


LOO 
.99901 
.98328 
.90608 
.86634 
.81394 
.74668 
.66738 
.644S9 
.40268 


1.00 

.99918 
.0S&7& 
.91733 
.88113 
.83248 
. 78783 
.68273 
.67163 
.42710 


1.00 
.99879 
.97977 
.88113 
.84726 
.78050 
.71822 
.62728 
.61421 
.37504 


1.00 
.9600 
.8400 
.6400 
.6775 
.6100 
.4376 
.3600 
.2776 
.1900 


.95 
1.00 
1.05 
1.10 
L20 
l.SO 
1.50 
2.00 


.22365 
0 

-.27762 
—.62060 
-L 66960 
-Z 95916 
-& 01227 
-56.6884 


.21041 
0 

-.30370 
-.70423 
-L857U 
-8.73944 
-11.8697 
-163.79 


.20534 
0 

-.21667 
-.64102 
-1.30I5S 
-2.S4882 
-5.64468 
-26.9081 


.0975 
0 

-.1026 
-.2100 
-.4400 
-.6000 
-1.2500 
-3.0000 




'7-l.-F=(l-A/^'** 
1 F=\—M* 



